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Converged approximate density functional calculations usually do not bind anions, 
due to large self-interaction error. But Hartree-Fock calculations have no such prob- 
lem, producing negative HOMO energies. A recently proposed scheme for calculating 
DFT energies on HP densities is shown to work very well for molecules, better than 
the common practice of restricting the basis set, except for cases like CN, where the 
HF density is too inaccurate due to spin contamination. 
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Anions and radicals are important for many applications, including environmental 
chemistry, semiconductors,^'^ fuUerene chemistry, ^"^"^ charge transfer, and solar cells. ^^'^^ 
Recently, electron affinities of biological species become of great interest, especially in studies 
of low-energy electron DNA damage. Low-energy electrons cause single-strand breaking, 
double-strand breaking, and supercoil loss in DNA even below the DNA ionization poten- 
tial. The electron affinity of DNA bases and base-pairs is important in determining damage 
mechanism. 

Density functional theory (DFT) has become a standard method for electronic struc- 
ture calculations in chemistry, and the standard functionals can be applied with standard 
basis sets to calculate electron affinities. The results are excellent, with mean absolute er- 
rors (MAE) below about 0.2 eV.^^ However, there is a theoretical fly in the computational 
ointment: Inspection of the orbital energies show that the HOMO of the anion is usually 
positive. This implies that, in principle, the calculation is unconverged.^^ If a sufficiently 
large basis set had been used, a fraction of the additional electron would ionize^^ and the 
HOMO drop to zero. This is due to the self-interaction error that all the standard density 
functional approximations suffer from. This error is especially large for anions, because of 
their additional electron. Self-interaction error produces an exchange-correlation potential 
that incorrectly decays exponentially in the asymptotic region, instead of decaying as — 1/r. 
For atomic anions, a large positive barrier appears in the Kohn-Sham potential (See Fig. 
1 of Ref 24) resulting in positive HOMO resonances. These metastable states are occupied 
and artificially bound by moderate basis sets (MBS), and so produce a positive HOMO. 
Because the positive barrer is often very wide, the total energy appears converged unless 
extreme basis functions are used. 

There have been strong discussions about this issue. ^^'^^"^^ Users find reasonable results 
with MBS for most cases, and ignore the postive HOMO. Purists regard all such calculations 
as unconverged, and so their results are suspect. Pragmatists will report results with the 
standard methods, but attach a caveat emptor footnote. The paradox has recently been 
addressed in several papers, ^'^'^^ which explain how accurate results can come from such 
unconverged calculations, but also suggest an alternative procedure that avoids the dilemma: 
Evaluate the density functional total energies on Hartree-Fock (HF) densities. We refer to 
this method as HF-DFT. 

Because HF is exact for one-electron systems, it has no self-interaction error, and its 
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HOMOs are bound, even for anions. Electron affinity calculations for atoms and their 
anions show excellent results with either method, ^^'^^ with MAEs about 0.1 eV, about half 
of that for ionization potentials. 
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FIG. 1. Calculated electron affinities of molecules in the G2-1 set (excluding CN) and atoms of the 
first two rows of periodic table plotted against experimental electron affinity. B3LYP energies were 
evaluated on self-consistent densities, while PBE energies were evaluated on HF densities within 
the AVTZ basis set. 

In the present work, we test the new procedure for the adiabatic electronic affinity of 
the molecules in the G2-1 data set.^^ All calculations are performed with TURBOMOLE 
6.2.29 The functionals used in DFT calculations are PBE,^^ B3LYP,3i-33 and PBEO.^^ We 
use Dunning's augmented correlation-consistent pVTZ (AVTZ) basis set.'^^'"^^ Structures of 
neutral molecules and molecular ions are optimized with self-consistent Kohn Sham-DFT 
and non-scaled zero-point vibrational correction with the same functional is added. For 
HF-DFT calculations, unrestricted HF calculations are performed on the DFT optimized 
structures. Based on these HF orbitals, the same functional used in structure optimization 
was selected for energy evaluation. We distinguish such calculations by HF-XC, where 
XC indicates the exchange- correlation approximation used. For all our calculations, SCF 
convergence was achieved with 10~^ eV or lower energy difference and with 10^^ root-mean- 
square density matrix element deviation. We exclude CN because the unrestricted HF 
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solution of neutral CN is strongly spin contaminated.^ This results in both a large energy 
destabilization and also poor HF densities. ^^'^^ 

In Fig. 1, we plot calculated versus experimental electron affinities, showing just how 
good the overall agreement is. In Table I, we report results for both methods for all molecules 
in the G2-1 set, but with averages excluding CN. Averages are reported in MAE and mean 
of errors (ME). We find once again excellent results for almost all molecules, with either 
conventional DFT with MBS or the HF-DFT method, with MAEs again about 0.1 eV. We 
find HF-PBE yields the best results overall, but differences are slight and not significant. 
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FIG. 2. Comparison of errors in electron affinities (eV) of molecules in the G2-1 molecule set 
excluding CN and atoms of the first two rows of periodic table. The colored dotted lines indicate 
the mean error of each method. All calculations use AVTZ basis set. 



In Fig. 2, we plot errors versus electron affinities, including also the results from atoms 
of Refs 24 and 27. We have excluded CN, where spin contamination of the HF wavefunction 
makes the HF density of the neutral too poor for accurate calculation. Spin contamination 
in HF for CN is > 50%, whereas no other molecule in the set has greater than 10%. It is 
suggested that calculations with spin contamination higher than 10% are not reliable.^" No- 
tice that switching the density from self-consistent to HF always either reduces the electron 
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affinity (sometimes increasing the error), or increases it by no more than 0.02 eV. In the 
case of PBE, the self-consistent electron affinities are mostly too large, a systematic error 
inherited from the local density approximation (LDA, sometimes called VWN'^^). This is 
reflected in the fact that the ME ~ MAE on the scale of the MAE in Table I. The reduction 
in electron affinities on using the HP density, which leads to subsequent reduction in MAE 
and large reduction in ME, shows that this is largely a self-interaction error in the density, 
not the energy functional. On the other hand, B3LYP is a hybrid functional with empirical 
parameters. It already cancels some self-interaction error, and has smaller MAE. But the 
ME is much smaller than the MAE, showing that its errors have random signs, i.e., much 
less systematic than those of PBE. Inserting the HP density overcounts the self-interaction, 
does not improve MAE, and even worsens ME. To check our interpretation of the effect 
of the hybrid, we applied another functional, the non-empirical hybrid PBEO'^^, and found 
results with the same trends, but higher MAEs (1.4 eV with or without HP densities). We 
also show the HOMO energies for the anions, in both HP and self-consistent calculations. 
All species except CI2 have positive HOMO in the DPT calculations, indicating their un- 
bound nature. A sufficiently large basis^® would reduce this value, but this effect may not 
be noticeable with any standard basis set. On the other hand, although the HP HOMO's 
are negative, they are not an accurate guide to the true electron affinities. Using Koopman's 
theorem, one may estimate electron affinity from the HOMO energy of anion, but both re- 
laxation and correlation effects are so large that Koopman's theorem is unhelpful here. Since 
the HP density is calculated upon geometries optimized from DPT, the resulting HP-DPT 
energy will not typically be a minimum in the HP-DPT potential energy surface. This shows 
further improvement may be made in HP-DPT by development of potential energy surface 
scan and optimization techniques. 

A useful tool for understanding these effects is the electron affinity density: 

riEAir) = n„i(r) - no{r) (1) 

where rio(r) is the charge density of the neutral, and n_i(r) is that of the anion. Pig. 6 
of Ref 27 shows this for the CI atom and anion, and how the HP electron affinity density 
is more compressed than that of using MBS and standard functionals. We plot the cross 
sections of electron affinity densities of NH along the molecular axis for different methods. 
In Pig. 3, the electron affinity densities of self- consistent and HP calculations are plotted 

5 



C3 



-1.5 




- (a.u.) 



. 


■ 

1 


0., 
0.08 




1 


0.06 




1 


0.04 






0,02 











1 


-0.02 






-0.04 




1 


-0.06 


1 




-0.08 
-0.1 



FIG. 3. Cross sections of electron affinity densities (anion - neutral) along the molecular axis of 
(a) PBE and (b) HF plotted for NH. 



for NH. NH, which has the largest deviation between the electron affinity error of PBE 
and HF-PBE, the electron affinity density of the two is nearly identical. Nonetheless, the 
self-consistent density is more diffuse than the HF density due to the electron leakage in 
the anion as shown in Fig. 4. To confirm this is not an artificial effect from the geometry 
difference in neutral and anion, we present the local vertical electron detachment densities, 
which are electron affinity densities using the same geometry for the anion and neutral. 

In summary, we have shown that the methods proposed in Refs 24 and 27 work just as 
well for the small molecules of the G2-1 data set, so long as the HF densities are close to the 
true densities. The results are equally good with HF-PBE as with B3LYP, perhaps slightly 
better. On the other hand, we found no case where the limited basis set approach fails. 
Presumably, the unphysical barrier holding the additional electron in for atomic anions^^ 
is sufficiently large for our molecules that standard basis sets show no sign (other than a 
positive HOMO) that the state being calculated is a resonance, rather than an eigenstate. 
We recommend that the HF-DFT method be applied more broadly for electron affinity 
calculations, especially for cases where DFT with MBS is believed to be inaccurate. 
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FIG. 4. Local electron affinity densities (anion - neutral) along the molecular axis (z-axis) in Fig. 
3. N is positioned at z = —0.99 and H is at z = 0.99. Regions in (a) are magnified into (b) and 
(c) for clarity. 
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TABLE I. Electron affinities (EA) of molecules and HOMO eigen values of anions in the G2-1 
set excluding CN (eV). All calculations were with the AVTZ basis set and using DFT optimized 
geometries. HE HOMO eigenvalues were evaluated from B3LYP geometries. 



A/Tnl 
iVlUi. 


EA 
expt. 


AEA(MBS) 
PBE B3LYP 


AEA(HF-DFT) 
PBE B3LYP 


£homo 
B3LYP HE 


CH 


1.24 


0.29 


0.04 


0.22 


-0.03 


-1.3 


2.1 


CH2 


0.65 


0.13 


0.06 


0.05 


0.02 


-1.3 


1.3 


CH3 


0.08 


0.00 


-0.05 


-0.04 


-0.09 


-1.7 


0.6 


NH 


0.38 


0.17 


0.02 


0.04 


-0.07 


-2.1 


0.1 


NH2 


0.77 


0.06 


-0.11 


-0.01 


-0.15 


-1.7 


1.3 


OH 


1.83 


0.12 


-0.09 


0.00 


-0.16 


-1.1 


3.0 


SiH 


1.28 


0.12 


-0.08 


0.10 


-0.12 


-0.8 


1.5 


SiH2 


1.12 


0.17 


-0.01 


0.13 


-0.06 


-1.0 


1.3 


SiHg 


1.41 


0.01 


-0.06 


0.01 


-0.05 


-0.3 


1.8 


PH 


1.03 


0.04 


-0.01 


0.02 


-0.01 


-1.1 


0.9 


PH2 


1.27 


-0.02 


-0.09 


-0.01 


-0.08 


-1.0 


1.2 


HS 


2.36 


-0.03 


-0.11 


-0.04 


-0.11 


-0.2 


2.6 


O2 


0.44 


0.00 


0.03 


-0.08 


-0.02 


-2.2 


2.4 


NO 


0.02 


0.27 


0.26 


0.16 


0.14 


-2.3 


2.5 


PO 


1.09 


0.18 


0.11 


0.14 


0.07 


-1.1 


2.0 


S2 


1.66 


-0.07 


-0.04 


-0.04 


-0.02 


-0.5 


2.2 


CI2 


2.39 


0.27 


0.40 


0.26 


0.38 


1.9 


4.7 


MAE 
ME 


0.00 
0.00 


0.11 
0.10 


0.09 
0.02 


0.08 
0.05 


0.09 
-0.02 


2.2 
-2.2 


* 

0.8 

* 

0.8 



Based on Koopman's theorem, mean absolute errors 
and mean errors are obtained by comparing the dif- 
ferences between — enoMo EA. 
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